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Abstract 

We study topologically stable non-Abelian global vortices in the U {N) linear sigma 
model. The profile functions of the solutions are numerically obtained. We investigate the 
behaviour of vortices in two limits in which masses of traceless or trace parts of massive 
bosons are much larger than the others. In the limit that the traceless parts are much 
heavier, we find a somewhat bizarre vortex solution carrying a non-integer U (1) winding 
number l/\/iV which is irrational in general. 
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1 Introduction 



Superfluid vortices appear in various condensed matter systems such as helium superfluid. They 
are global vortices in relativistic field theories [1]. When a global f/(l) symmetry is spontaneously 
broken, for instance by the order parameter (0) = v with a complex scalar field in the Goldstone 
model, a f/(l) Nambu- Goldstone boson appears. Then in general there appears a global string 
asymptotically ~ f e*^ winding around the vacuum manifold (or the order parameter space) 
U{\). The energy of global strings is logarithmically divergent in the infinitely large space, 
unlike local vortices with gauged U{1). So less attention have been paid to global vortices as a 
candidate of cosmic strings compared with local counter part. When the net number of strings 
is zero, for instance, a pair of string and anti-string or string loops, their energy is finite and can 
be considered as cosmic strings [2ll3]. Axion strings are such objects. 

Natural non-Abelian extension of U{1) vortices is U{N) vortices in the U{N) linear sigma 
model for which the order parameter is extended to an by complex matrix ($) = wIatIi] 
Such non-Abelian vortices are expected to form during the chiral phase transition in QCD, in 
which case the field $ of a 3 by 3 matrix (A^ = 3) is a condensate of quark-anti-quark ($) ~ (qq)- 
At the chiral phase transition [l], the chiral symmetry SU{N)-[^ x SU{N)Yi is spontaneously 
broken down to its diagonal symmetry SU{N)y. According to this breaking, massless Nambu- 
Goldstone bosons appear as pions (or more generally mesons). At the same time, the axial 
symmetry U{1)a is also spontaneously broken and the rj' meson appears. However U{1)a is 
explicitly broken by the axial anomaly at zero temperature, giving a mass to the r]' meson. It 
has been argued that the axial anomaly might disappear and U{1)a is approximately recovered 
at high temperature [5]. Although there is still an ambiguity if it occurs below the temperature 
of the chiral phase transition, let us consider such a situation. Then the breaking pattern is 
U{1)a X SU{N)i^ X SU{N)ji SU{N)y apart from the discrete symmetry, and the vacuum 
manifold is 

'-Em^^^u(Nu. (1.1, 

The first homotopy group 7ri[U{N)A] — Z is non-trivial and so there exist topologically stable 

^ Other non-Abclian global vortices appear in the B-phase of ■^He in which symmetry is broken as SO{3)s x 
50(3)l X C/(1) ^ 50(3)s+L- In this case the corresponding vacuum manifold is C/(l) x S0{3) and the first 
homotopy group is 7ri[;7(l) x SO {3)] ~ Z © Zj. 
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vortex-strings in this breaking. The simplest vortex appearing is a U{1) vortex-string called 
the rj' string asymptotically given by $ ~ ve^^lj\f, which winds around U{1)a once (27r) [6l 
[Tjl^ However this is not the minimum vortex-string because there exists a smaller loop in the 
vacuum manifold (11.11) : The minimum string is a non-Abelian string asymptotically given by 
$ ~ t'diag (e*^, 1, ■ ■ ■ , 1) = w e*^diag ^e*"*^^, e~*^, ■ ■ ■ , e~*^ j , which winds U{1)a as well as the 
SU{N)a [9]. The important is that this string winds 271 /N of U{1)a and therefore its tension 
is 1/A^ of that of the U{1) T^'-string. Such vortices with fractional f/(l) winding number often 
appear in various condensed matter systems such as Bose-Einstein condensates and certain types 
of superconductors, and are called "fractional vortices" If] 

At low temperature the axial anomaly induces the periodic potential in U{1)a- The U{1)a is 
broken to Zat and there appear disconnected vacua, in each of which the rj' meson gets mass. 
A U{1) ?7'-string is accompanied with A^ domain walls and the total configuration becomes A^ 
domain wall junction with a string at the junction line [10]. Balachandran et. al have discussed 
a possible role of such an object in the early universe [TO] . 

The presence of a non-Abelian vortex breaks the SU (A^)v symmetry of vacua to its subgroup 

SU{N — l)v X f/(l)v and consequently there appear further Nambu-Goldstone modes CP^^^ ~ 

SU{N)/[SU{N — 1) X f/(l)] which are orientations in the internal space [TT]. This idea was 

brought from the local U{N) vortices [12] for which U{1)a and SU{N)i^ are gauged. However 

there exists a crucial difference between global and local U (N) vortices. The Nambu-Goldstone 

modes CP^~^ of local U{N) vortices are localized around the vortex and become the moduli (or 

collective coordinates) of the vortex [12] while those of global U{N) vortices are not localized but 

spread to infinity (or the boundary of a finite space). Having this in mind, an inter-string force 

between two parallel global U{N) vortex-strings with different orientations in the internal space 

has been calculated recently [131 E]. The force depends on the relative orientation: it reaches the 

^ Brandenberger et. al have discussed non-topological strings, called the pion strings, using massless Nambu- 
Goldstone (pions) SU{2)a [6l[7l[8]. They are topologically unstable because of tti[SU{N)] = 0. In this sense 

non-Abelian strings below are made of both rj' mesons and pions. 

For instance in the polar phase of a spin 1 spinor Bose-Einstein condensate, [/(!)$ x S0{3)s is spontaneously 

broken down to C/(l)$+5 x: (Z2)$+s with x denoting a semi-direct product. Then the vacuum manifold is 

M ~ [[/(!)$ X SO{3)s]/[U{l)^+s X (Z2)<i,+s] - [U{1) x S^]/Z2 |15J and the first homotopy group tti{M) ~ Z 

supports half quantized vortices [16| . Similarly to this, 1/3 quantized vortices exist in the cyclic phase of a spin 

2 spinor Bose-Einstein condensate [T7] . 
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maximum when two strings wind around the same component of the vacuum expectation value 
[for instance diag (e*^, 1, ■ ■ ■ ,1) and diag (e*^, 1, ■ ■ ■ ,1)], and it vanishes when the two strings wind 
around different components of the vacuum expectation value [for instance diag (e*^, 1, 1, ■ ■ ■ ,1) 
and diag (1, e*^, 1, ■ ■ ■ ,1)]. This result implies that a U (1) r^'-string is marginally decomposed into 
N pieces of non-Abelian strings: e^^l^ diag (e*^, 1, 1, ■ ■ ■ , l) + diag (1, e*^, 1, ■ ■ ■ , !) + ■ ■ ■ . Such 
decomposition necessary occurs at finite temperature where the free energy is minimized instead 
of the energy. Therefore at low temperature with the axial anomaly, a domain wall junction [TU] 
is unstable, because a string is pulled by each of domain walls and is decomposed into 
non-Abelian strings to each of which one domain wall is attached. In the end each piece is 
pulled to infinity in each of A^ directions^ 

In this paper we study the purely string solution without domain walls in the linear sigma 
model without the axial anomaly. Numerical solutions themselves were previously obtained in 
Here we study profile functions of solutions in much more detail with more accuracy. By 
using the relaxation method, we numerically determine the shooting parameters of the solutions 
up to fifth order. We investigate the dependence of the profiles of the vortex to the parameters 
in the linear sigma model. We also study the two limits in which the masses of traceless or 
trace parts of the massive bosons in the linear sigma model are much larger than the others. In 
these limits, the model reduces some nonlinear sigma models. In the limit that the trace parts 
are much heavier, the equations for the profiles of the U{N) vortex become sine-Gordon-like. 
The solution remains to be regular. On the other hand, in the limit that the traceless parts 
are much heavier, we find somewhat surprising solution; the U{N) vortex solution reduces to 
a singular U{1) vortex with the U{1) winding number 1/Vn which is irrational specifically for 
A^ = 2, 3, 5, 6, ■ ■ ■ . In general, profiles of U{1) vortices with a non-integer U{1) winding number 
(< 1) are of course singular. Interesting is that such an irrational vortes^^ naturally appears in 
a particular limit of a regular non-Abelian vortex solution. As far as we know such a vortex has 
not been reported yet in the literature. 

This paper is organized as follows. In Sec. [2] we review the U{1) global vortex solution in the 
Goldstone model. We give numerical solutions and study their asymptotic behaviours. In Sec. [3] 
we study the non-Abelian vortex solution in the U{N) linear sigma model. We derive profile 
^Stable domain wall junctions or networks exist in gauged U(N) linear sigma models with appropriate masses 
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functions of the minimum U{N) vortex with the U{1) winding number for = 2, 3, ■ ■ ■ , 10 
and determine the shooting parameters up to fifth order. We also discuss various hmits by 
sending some of masses to infinity. We find that in a particular limit the non-Abelian U{N) 
vortex reduces to an Abelian vortex with the irrational U{1) winding number l/\fN. Sec. H] 
is devoted to conclusion and discussion. Behaviors of vortex solutions in the large limit is 
discussed in Appendix. 

2 Global U(l) Vortices 

2.1 The Goldstone Model and Vortex Solutions 

Let us begin with giving a review on a global vortex-string solution in the Goldstone model with 
a complex scalar field 

/: = |9,0p-A(|0p-^2)'. (2.1) 

We choose A > and f ^ > for stable vacua with broken U{\) symmetry. The scalar potential 
is like a wine bottle, so we have an vacuum space with radius |0| = v. When we choose a 
vacuum cj) = v and consider small fluctuations as (j) = v + "^^^^ {¥^1,2 ^ the Lagrangian in 
terms of the small fluctuations up to quadratic terms is of the form 

^^'^ = lid.^if + 1(9,^2)' - 2XvVi- (2.2) 
This shows that ipi is massive and is a massless Nambu-Goldstone mode: 

mj = 4A^;^ mj = 0. (2.3) 

The equation of motion of reads 

9^a'^0 + 2A0(|0|'-t;') = 0. (2.4) 

A global vortex-string extending linearly to the xs-axis is obtained by solving (12.41) with an 
axisymmetric vortex ansatz in the cylindrical coordinates xi + ix2 = tc*^, given by 

(j){r,e) =ve'''^f{r), k e Z, (2.5) 
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with the boundary conditions 

hm /(r) = 1, hm/(r) = 0. (2.6) 

r— >oo r— >0 

Plugging the ansatz f l2.5p into Eq. (12.41) . we get a second order differential equation 

/" + 7-^-^/(/'-i)=0- (2-7) 

Numerical solutions for k = 1,2, ■ ■ ■ ,10 with the boundary conditions (12.61) are plotted in the 
left panel of Fig. [H 




Fig. 1: Left panel shows the profile functions and the corresponding energy densities S{k) (solid 
lines) and Spot{k) (broken lines) are plotted in the right panel for k = 1,2, ■ ■ ■ ,10 and mi = 1. 



The energy of the vortex solution can be expressed as 

r + 



E = 271 I drrS = 2txv'^ / drr 
'o io 



(2. 



where the first two terms come from the derivative of the field and the last term is from the 
scalar potential. Since / — ^> 1 as r — oo, the kinetic energy logarithmically diverges. So the 
energy of the vortex-string consists of a finite part and a logarithmically divergent part as 



E{k) = Ederik) + E^,{k) + Ep,,{k) 



where we have defined 

Ederik) 

Edivik) 



poo 

2ttv^ / drrf'^, 
Jo 

r"00 



271 / drrS, 



2 roo 

Pot = ^ drr{f-lf 



2 

oo _f 2 



9 / L 

2Txv'^k'^ I dr — = const. + 2'Kv'^k'^ lim log — , 

r L^oo To 



(2.9) 

(2.10) 
(2.11) 
(2.12) 



Here we have introduced an IR cut-off L which is the size of the system and ro ~ m]^ is a typical 
size of the U{1) global vortex such as | / — 1| ^ 1 as r ^ tq. We can analytically calculate Ep^t 
as follows. Let us first introduce a dimensionless coordinate p = rriir. The equation becomes 
independent of the coupling constants 

F[f{py,k] ^f"+^--^~\f{P-l)= 0. (2.13) 
p p z, 

Then we see that -Epot is independent of the scalar mass mi. 

EMk) = ^ f drrif - if = 2nv' dpp^-^^^. (2.14) 
To calculate this, let us make a trick 

+ f(/'-l)'- (2.15) 
By using this equation, we can bring Epot in the following form 

oo 

= nv^k'^. (2.16) 



This formula is called the Derrick-Pohozaev identity in the literature [TH] . To derive the rightmost 
equality, we have used the asymptotic behaviours / oc p'^ at p ^ 1 and / = 1 — a2/p^ + 0{p~^) 
for p ^ 1, which will be obtained in the next subsection. We plot the total and potential energy 
densities S and Spot, respectively of our numerical solutions in the right panel of Fig. [H We also 
have numerically checked Epo^/ {jiv'^k'^) = 0.999 + (9(10~^). We can show that E^^r is also finite. 

Note that higher winding solutions {k >2), especially co-axial vortices as we assumed above, 
are not stable. Since the energy of the co-axial vortices is almost proportional to k^, distant 
vortices are energetically preferred (/c^ > k for k > 2). Therefore, the above static solutions with 
k > 1 are artifacts of our co-axial ansatz and boundary conditions. 

2.2 Asymptotic behaviors at r ^ and r ^ oo 

Let us investigate asymptotic behavior at mir = p ^ 1 where the profile function is very small 
I/I ^ 1. We expand the profile function as 

oo 

/(p) = ^a„p^ (2.17) 

n=l 



= 2p'f'F[f;k] 



dp 



p'r-k'p-p-u'-if 



Epotik) 



TTV 



-f 



12 



if 



2 r-2 



k'f 
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By substituting this in Eq. (12. 131) . we find that the first non-zero coefficient is that of 

afc = hm^ (2.18) 

which is sometimes called the shooting parameter and which may be determined by making use 
of numerical solutions. Such parameters are important since we can uniquely determine solutions 
with the parameters. In the minimal winding {k = 1) vortex, we got ai = 0.4123772 and 

fk=i(p) = aip -p^ + — - — + C p^ . 2.19 

In principle, one can infinitely increase accuracy of the approximation with the unique parameter 
ai. Although we assumed p <^ 1 in the beginning, we can reach at po > 1 with a good accuracy 
by increasing the order of expansion. 

Let us next discuss asymptotics at r — >■ 00. We expand solution in the following way 

/(P) = l-E4' (2-20) 

Up 

hi being constant. Unlike the expansion parameters {aj} in Eq. (12.171) . the expansion parameters 
{hi} can be precisely determined. To this end, we insert fl2.20p into F[f\ k] given in Eq. fl2.13l) 
and determine hi by comparing terms order by order. Then we get 6odd = and 

62 = 64 = -k'^ (8 + , fee = -^^ (l28 + 32^;^ + fc^) , ■ ■ • (2.21) 
2 2 

We can go on up to order which we desire. 



3 Global U{N) Vortices 

Let us study non-Abelian global vortices. Basic analysis on them have been done in [11]. Here we 
are going to push forward the analysis of [11] on vortex solutions by investigating the equations 
of motion in much detail. 



3.1 The U{N) Linear Sigma Model and Vortex Solutions 



The model that we consider here is a natural extension of the Goldstone model given in Eq. (12. ip . 
It is the SU{N)i^ x SU{N)^ x U{1)a linear sigma model for an x complex matrix of scalar 
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fields given by 

C = Tt [9^$t5M$ _ A2($t$)2 + ^2$t$] _ X, (Tr[$t$])2 - (3.1) 

Al + A2) 

where the last constant term is introduced for the vacuum energy to vanish. For a stability of 
vacua, we consider the parameter region /i^ > 0, A2 > and NXi + A2 > 0. The chiral symmetry 
SU{N)i^ X SU{N)^ and the axial symmetry U{1)a act on $ as 

$ ^ e''gL<^gR, (e^^ Ql, Qr) e U{1)a x SU{N)^ x SU{N)^. (3.2) 

However, unlike the usual case in the absence of U (1)a broken by the axial anomaly, the structure 

of discrete symmetries becomes somewhat complicated in the presence of U{1)a- Here we explain 

it in detail. First the group action G on $ is not Eq. (13. 2p itself but is given by 

^_ Uil)AxSUiN)^xSUiN)n ^ ^gg^ 
X 

where the following Zat x Z^r action does not act on $ and therefore is removed from G: 

{uj-''-\u;H^,uj'1n) G U{1)a x SU{N)i^ x SU{N)r, 

^^^2m/N^ (A;,/ = 0,1,2,--- ,iV- 1). (3.4) 

For later use let us redefine these discrete groups as Z^ x Zjy — {'^n)v x {'Zn)a with 

(Za)v : (1, uj'^In, uj~'1m) G U{1)a x 5f/(iV)L x SU{N)j, (3.5) 
(Za)a : (a;-2^^'=l;v,^'l7v) G U{1)a x 5f/(iV)L x 5t/(iV)R. (3.6) 

By using the symmetry G, any vacuum can be transformed into the form 

<*> = "^- "' = m^y 

We can consider this vacuum without loss of generality. All other degenerate vacua are obtained 
from this by the G transformations. The symmetry G is spontaneously broken down to the 
isotropy group 

^ ^ SUiN)y X (Za)a ^ SUjNh ^ 

Za X Ztv (^Af)v 

where Z^r x Z^ is the one of Eq. (13.41) which can be rewritten by Eqs. (13. 5p and (13. 6p . and 
SU {N)y is given by 

SUiN)y : (1, g, g^) G U{1)a x ^f/(iV)L x SUiN)^. (3.9) 



When the vacuum ($) is transformed by G, the isotropy group H is also transformed by a 
similarity transformation, but all of them are isomorphic to the original isotropy group fl3.8p . 
Therefore the vacuum manifold (the order parameter space) can be written as a coset space: 

G ^ (^(1)a X SU{N)^ X SU{N)j,)/{'Ln x Z^y) 

H SU{N)y/{ZN)y 

which is eventually a group manifold (because if is a normal subgroup of G). Since the first 
homotopy group of the vacuum manifold 

7ri[f/(iV)A] ~Z (3.11) 

is non-trivial, it admits topological vortex-string solutions. 

In order to find the mass spectrum, let us perturb $ with small fluctuations as 

$(x) =t;l,v + 0(x)l7v + Xa(x)T", ia = l,2,--- -1), (3.12) 

with the generators T'^ of SU{N) (Tr [T'^T*] = 5'^^). Then it is turned out that imaginary parts 
of and Xa are massless Nambu-Goldstone modes parametrizing the vacuum manifold U{N)a 
in Eq. (13.101) . while their real parts are massive: 

mJ = 2/.^ ml = ^^^^^^^^ = AX,v\ (3.13) 

The original coupling constants are expressed by these masses as 

2 2 2 2 

The scalar potential can be rewritten by the dimensionfuU parameters v, and as 

^ = (Tr [$^$ - ^'^1^])' + §Tr [($t$)^j , (3.15) 



where we have introduced a notation (X) = X — -^^^l^r for X x X matrix X (note (Iat) = 0). 



Let us introduce a ratio of the two masses 



We will see that r determines all properties of non-Abelian global vortices. 



'Instead of this, k = A1/A2 was used to parametrize solutions in [TT] . 



Let us construct vortex-string solutions in this model. Firstly, one may consider the following 
simple boundary condition 

lim $(r,^) = t;e^^l;v (3.17) 

A natural ansatz with this boundary condition is 

^r,e)=ve''f{r)lN, (3.18) 

with boundary conditions /(O) = and /(oo) = 1. The equation of motion for the profile 
function /(r) is identical to that of the global U{1) vortex with k = 1 

r + 7-^-^/(/^-l)=0- (3-19) 

So the solution is just embedding of the global [/(l)-vortex of Eq. (12.131) . It is important to 
observe that neither nor appears. Furthermore, only the overall U{1) phase winds once 
when we go around the vortex solution. This Abelian vortex solution is called the t]' string [7]. 
However, as we will see, this solution is not minimal [11] and is broken into A^ minimal solutions 

m- 

A minimal winding vortex configuration can be obtained by taking the following boundary 
condition 

lim <l>(r, 9)=v diag(e*^ 1, ■ ■ ■ , 1) = w e^^diag f e^^^ e"^^, ■ ■ ■ , e"^^) . (3.20) 

r— >oo V / 

The key point here is that the overall f/(l)-phase winds only 211 /N around the vortex and the 
rests 27t{N — 1)/N or —2'k/N are inside the non-Abelian group SU{N)x- As a consequence, the 
tension of the vortex is 1/A^ of the above Abelian solution fl3.18p . Because the vacuum space 
is U{N)x, the above vortex is called a non-Abelian global vortex or more specifically a global 
U{N) vortex. Corresponding vortex ansatz is of the form 

$(r, 9)=v diag {e''f{r),g{r), ■ ■ ■ , g{r)) , (3.21) 

with the boundary conditions 

lim / = lim ^ = 1, lim / = 0, lim g' = 0. (3.22) 

r— >oo r— >oo r— >0 r~*0 

Since / is the profile function of a winding scalar field, it must go to zero at the origin. On the 
other hand, g does not have to vanish. 
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For a fixed 9 = Oq, the field $ behaves at the boundary as = 9o,r ^ oo) = v diag(e*^°, 1, ■ ■ ■ 
where the isotropy group is not the same with the one in Eq. fl3.8p . Instead, it is obtained 
from H in Eq. ([HS]) by a G transformation (e^^,^o,^o) e U{1)a x SU{N)i^ x SU{N)p^ with 

= diag ^e^'w-^o, e'^aw-, ■ ■ ■ ,e~*2w-j, and therefore it is isomorphic to if = SU{N)y / {'^n)y 
in Eq. (13.81) . Around the vortex where g{r) differs from /(r) (especially at the center $(0) = 
f diag(0, (7(0), ■ ■ ■ , (7(0))), the isotropy group Hq^q^ = 5't/(A^)v,0o/(^Af)v,0o further broken to its 
subgroup [SU{N — l)y x [/(l)v]/(ZAr)v,eo- This gives rises to further Nambu-Goldstone modes 

m 

In the case of the local U{N) vortices for which U{1) and SU{N)i^ are gauged [12], the corre- 
sponding modes are called the orientational zero modes of a vortex. However, these are non- 
normalizable for a global U{N) vortex (in the infinite space), unlike the local U{N) vortex, 
because the isotropy groups Hq depends on ^ = ^0 ^-^d differ from each other, and consequently 
the wave functions of the Nambu-Goldstone modes CP^^^ spread to infinity. 

Let us investigate concrete profile functions of vortex solutions. The Hamiltonian density and 
the energy densities in terms of the profile functions /(r) and g{r) is given by 

n = 2Wr8, (3.24) 



£ = f + L + {N-l)g'^ + V, (3.25) 
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^ = '§U'+iN-l),^-Nf+''-^i^(f-,^f. (3.26) 
By minimizing the Hamiltonian 7i, we get the equations of motion for / and g 

r + 7 - ^ - [f + 1)/ -N)- ^^^^^/ if - = 0, (3.27) 

^" + 7 - if + (iV - 1)/ -N) + ^9 if - = 0. (3.28) 

The third term in the left hand side of Eq. (13.271) is typical for global vortices. This leads to 
logarithmic energy divergence. With respect to a dimensionless coordinate p = rrirpr, the above 
equations can be written in the following forms 

:FM[f,g,r] = f +---, — — = 0,(3.29) 

QN[f, 9:r]^ g" + '--^ if + (iV - 1)/ - iV) + (/^ - = 0. (3.30) 
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It is clear in this form that solutions depend on r only. We should solve these ordinary differential 
equations with boundary conditions fl3.22p with replacing r by p. 

Let us see the non-Abelian vortex in a special case r = 1. In this case, ^Ar[/, 1] = can be 
solved hy g = \ while the other equation J-'nU^ g = 1;t = \] = Q \s the same as fl2.13p . Therefore 
in the case of r = 1 the Abelian vortex is embedded into $ as a non-Abelian vortex solution. 
Qualitative behaviors of vortex profiles change at r = 1. As we take r to larger than 1, the 



f ,g B 




2 4 6 8 10 0123456 



Fig. 2: Configurations (solid lines for /(r) and broken lines for g{r)) in the left panel and energy 
densities in the right panel for logr = 0, 1/2, 1, 3/2, ■ ■ ■ , 9/2, 5 for N = 2. 



f'S s 




02468 10 02468 10 

Fig. 3: Configurations (solid lines for /(r) and broken lines for g{r)) in the left panel and energy 
densities in the right panel for — logr = 0, 1/10, 1/5, ■ ■ ■ , 9/10, 1 with N = 2. 

value g{0) gradually goes down toward zero and the energy density becomes sharp (but remains 
regular for finite r), see Fig. [21 On the other hand, g{0) grows and the vortex remains regular 
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and finite when we take r smaller than 1, see Fig. [3l As we will see below, g{0) is an important 
value. 

In order to see dependence of the non-Abelian global vortex, we plot two figures in Fig. H] 
by changing A^ form 2 to 10 with r being fixed to 1/2 and 2. 



f'S s 
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0123456 0123456 

Fig. 4: Configurations (solid lines for /(r) and broken lines for g{r)) in the left two panels and 
energy densities in the right two panels, with A^ = 2, 3, ■ ■ ■ ,10. The upper(lower) two figures are 
for T = l/2(r = 2). 

3.2 Asymptotic behaviors 

Let us investigate asymptotics of the non-Abelian global vortex. This has been studied in [TT] and 
here we want to extend it with a better accuracy and determine some fundamental parameters 
(shooting parameters) accompanied by the differential equations (13.291) and (13.301) . 

We start with study on asymptotics at r ^ max{mT^, m^^}. As usual, we expand the fields 
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oo oo 



/, g in the following 

/ = l-^a,r-\ g = l-J2b,r-\ (3.31) 

1=1 i=l 

By plugging these into Eqs. (13.291) and (13.301) and equating them with order by order, we can 
analytically determine the coefficients. It is immediately turned out that a2j+i and fe2j+i vanish. 
The leading order terms are of the form 

/.I 1 (i,,^^), ,.1 1 (i,_i,). (3.32) 

Note that / is always less than 1 while (7 ^ 1 for 1 ^ r. As observed in [11], the behavior of g{r) 
depends on the coupling constants Ai, A2. With respect to the physical masses and m<^, we 
now understand that the behavior depends on the ratio r of the two masses. The higher order 
terms are determined as 



(-1 + iV)(-l + m)ml + 2(-l + N)mlml + (1 + %N)m 



4 

04 = ^ — — , 4 — 4 — ' — ~ — ~ ^ ' — ^' (3.33) 



W)ml - 2m\m\ + (1 + 8iV)-^ 



X 



ae = 2Ar3m6m6 " ^)(^ + ^("^S + \^\N))m\ + (3 - (86 - %'iN)N)m\m\ 



(3.34) 



-(3 + (5 - m)N)m\m\ + (1 + 32iV(l + 4iV))mJ) , (3.35) 

2 2 

^^^|-^( (1 - 56A^ + 152Ar^) m\ + 2 (-1 + 12Ar + ^^N^) m\m\ 

+ (1 + 32A^ + 128A^2) m^) , (3.36) 



Let us next consider asymptotics at r ^ min{m^^, m^^}. We expand the fields by 

00 00 
] = Y,c,t\ g = Y,diT\ (3.37) 

i=0 i=0 

By plugging this into Eqs. (I3.29P and (13.301) . one can determine c^m = d2m+i = 0. The leading 
order of the approximation gives us 

f = c,r + 0{r^), g = do + Oir'). (3.38) 

The two parameters Ci = /'(O) and do = g{0) cannot be obtained analytically, but we can do it 
numerically. Note that they depend on r and A^, see Table [TJ 
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Table 1: Numerical data for the 
has analytic values, given in Eq. 
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Higher order terms are determined by Ci and do: 

N{l-d^o)ml+{ml-ml)dl 
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d4 
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8A^ 



■ dr, 



-Nml + (AT - 1) (mj - m^) dl 



16iV 



Cl, 



256iV2 



mcl ( 



2 2 \ 



+ (iV (3dg -I) ml- Ml {ml - mj)) (iV (rfg - l) mj - rfg (mj - m^)) 



Cl 



768iV2 



(iV2 _ 6(iV - l)iVd2 ^ _ ^)2^4^ ^ 

-2(iV - l)dl (-3iV + (-5 + 3iV)rf2j mjmj + (5 + (iV - 6+)iV)d^m: 
+lQNcl [ml + {N - l)ml) 



X 



(3.39) 
(3.40) 

(3.41) 



(3.42) 
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3.3 Heavy Particle Limits 



In this subsection, we are going to study how the vortices behave in two regions i) ^ 
and ii) ^ m^. Before doing it, remember that the hnear sigma model (13.11) becomes the 
chiral Lagrangian in the hmit m^, oo. This is because, after the heavy fields are integrated 

out, only massless NG modes ($ G U{N)) survive in a low energy theory. Our purpose of 
this subsection is to clarify effects of lightest massive fields to the chiral Lagrangian. For later 
convenience, let us rewrite the equations of motion (13.271) and (13.281) as follows 

f" f 1 ( n" n' \ m? 

^ + ^-- + (iV-l) - + - -^(f + (iV-lk'-iV) = 0, (3.43) 
f rf \g rgj 2 ^ 

f" f 1 n" n' 777^ 

T + 4 - 4 - - - - - ^ (/' - = 0- (3-44) 

J rj g rg Z 

A point is that appears only in the first equation while m^^ is in the second equation. 

Similar limits (with some particles being infinitely heavy) have been studied recently for a 
local U (N) vortex, for which U{1) and SU (A^)l are gauged, at the critical coupling (with a certain 
relation between gauge and scalar couplings) |20j . 

3.3.1 oo limit (r 0) 

Let us first consider the mass of the trace part of $ is much greater than that of traceless part, 

> m^. (3.45) 

In other words, we consider a limit r — ^ by sending — > oo with m,^ being fixeqj. Then we 
integrated out the heavy modes with mass m^. The first term in the scalar potential in Eq. (13.151) 
becomes very high and sharp, so that the scalar fields are squeezed into the following manifold 

Tr[$t$ _^;2i^] = 0. (3.46) 

To solve this condition, let us expand $ by 

$ = (^-f, f = (l^/yiV,T\T2,- ■■ ,T^'-^), (3.47) 



^ The same limit can be realized by taking with being fixed. However, the massless limit is tricky 

because the vortex infinitely spreads out and dilutes. 
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with A^^ complex vector (p = {ipo, ipi, ■ ■ ■ , ipf^2_i). Then the condition becomes 

|<^P = A^^;^ (3.48) 

Therefore the hnear sigma model (13.11) reduces to 5"^^^"^ non-linear sigma model. The scalar 
potential takes the form 

Vo = ^Tr [($t$)^] = gxr [($t$ _ ,^1^)^] . (3.49) 
The vacuum remains as U{N)p^. By using the chiral symmetry, we can choose 

= ^/Nv, ifa>l = ^ $ = vIn- (3.50) 

Let us next consider a vortex-string solution in this model by taking the same diagonal ansatz 
(13.211) as before. However, / and g are no longer independent because of the condition (I3.48p . 
They should satisfy 

f{rf + {N-l)g{rf = N. (3.51) 

With this at hand, we are aware of an interesting phenomenon that the shooting parameter 
g{0) = do approaches to an analytic value as 



lirndo = lim,(0) = y^. (3.52) 

Here we have used /(O) = 0. Our numerical result matches with this result, see Table [H The 
same thing can be found from the view point of equations of motion (13.431) . By dividing its both 
hands by and taking the limit ^ oo, we get Eq. (I3.5ip again. 

To simplify the other equation (13.441) . we may rewrite the fields by 



N 

/ = ViVcose, ^7=\/^^sine. (3.53) 



The reduced model is like the Sine-Gordon model: the potential is 

r^O: V = ^liE^4^. (3.54) 



The corresponding equation of motion is of the form 

^„ 6' 1 ^/2 . iV(iVcos2 - 1)\ 
r = 0: 0" + — + -sin2e^ + m2^— ^ = 0. 3.55 

r 4 ^ N - 1 J ^ 

We solve this with the boundary conditions 

lim6 = — , lim = arccos f ^= | . (3.56) 



We numerically solved this and the solution is shown in Fig. O 
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20 f,g 
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Fig. 5: Profiles of and corresponding f,g for N = 2 case (m^ = !)• 



3.3.2 — > oo limit (r oo) 

Next let us investigate the other limit 

:$> m^, (r oo). (3.57) 

We send infinity with being kept finite. Then A^^ — 1 real scalar fields in $ become 
infinitely heavy and we integrated them out from the theory. The scalar potential (13.151) gives 
us the following condition 

Tr [($"^$>^j =0 ^ (<l>"f$) = ^ <l>"^$ oc 1^. (3.58) 

We can solve this by 

$(x) = s{x) U{x) with see, U e SU{N)a. (3.59) 

This decomposition is up to the identification (s, U) ~ {uj'^s, u' -^U) with uj = e^^'/^ {k = 
0, 1, 2, ■ ■ ■ ,N — 1). Thus the theory in this limit is the SU{N) chiral Lagrangian coupled with a 
complex scalar field s. The kinetic term is given by 

i^™ = N\d^s\^ + \s\''Tr{d^Ud^U^) 

= N{d^af + (j^i:i{dJjd^U^) (3.60) 

with s = ae*" and U = e*"?7 G U{N) ^ [SU{N) x f/(l)]/Zjv. From the second expression, we 
see that the metric of the target space is a cone over U{N). The scalar potential becomes in this 
limit as 
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Thus the vacuum manifold remains to be U{N)a- The vacuum expectation value v of |s| gives 
the pion decay constant = 16|sp of the U{N) chiral Lagrangian. 

We are ready to consider vortex-string solutions in the limit. With respect to the profile 
functions f,g, the condition (13.581) forces us the following condition 

f{r)=g{r). (3.62) 

In terms of s and U, the non-Abelian vortex solution (13.211) can be expressed as 

s = e''^f{r), U = diag ^e'''^'^'"' , e"*^, ■ ■ ■ ,e"*^j . (3.63) 

Note that the overall U{1) winding of the non-Abelian vortex is 1/A^ as before. The condition 
(13.621) explains the behavior of g{0) which tends to go down toward zero as r is sent to oo, see 
Fig. [2] and Table [H Because of this behavior, the scalar fields $ vanish at the origin entirely so 
that the full chiral symmetry is recovered at the center of vortex. Furthermore with recalling the 
discussion around Eq. (I3.23p . we find that no second symmetry breaking occurs in the presence 
of the vortex solution (I3.63p . because the isotropy group is Hg even for finite 0) which is 
isomorphic to H = SU{N)y / in Eq. (13. 8p . Therefore, we find that the solution loses the 
internal orientations of CP^^^. Eq. (13.440 is automatically satisfied while Eq. (13.430 reduces to 

r = oo: r + f -^-^/(f-l)=0. (3.64) 

Let us compare this with Eq. (12. 7p for the Abelian vortex string. Interestingly, the non-Abelian 
global vortex in the — ^ oo limit can be identified with the Abelian global vortex with a 
non-integer U{1) winding number 




(3.65) 



which is smaller than unity, and can be an irrational numbers for A^'s which are not able to be 
expressed by squared integers. Of course, U{1) vortices with a non-integer U{1) winding number 
are generically singular because ~ e*^/v^ (at r ~ 0) is not single valued. Such singular 
solutions for F[f{p),^] = in Eq. (I2.13P with the boundary conditions (12. 6p are shown in 
Fig. El 

There is no smooth interpolation in the solutions between < oo and = oo in the 
following sense: Eqs. (I3.43P and (13.440 with finite always provide / ~ r to leading order in 
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Fig. 6: The left panel shows the profile functions and the corresponding energy densities S in 
Eq. ([2SD are plotted in the right panel for k = l/\fN with = 1,2,--- ,10. The Abelian 
solution of = 1 is only regular but the others are singular. 

r ~ (see Tabled]), and keep the energy density finite, while in the limit where = oo the 

1 2( — 1^ 

resultant Eq. ( 13.64^ gives f for r ~ and energy density gets divergent as r 7 at 

the center of vortex. Furthermore, the profile function g (in this limit g = f) does not satisfy 
the original boundary condition g'{0) = 0. Instead, it is replaced by g{0) = 0. 

Although this singular solution is an artifact appearing only in the limit where we have 
discarded the heavy modes completely, it reasonably accounts for the fact that the energy profile 
of the non- Abelian vortex becomes very sharp and finally looks singular when ^ m^, see 
Fig. O Taking into account such heavy modes, the singularity is smeared. 

4 Conclusion and Discussion 

In this article, we have investigated non- Abelian global vortices in SU{N)i^ x SU{N)-£i x U{1)a 
linear sigma model in detail. We push forward the analysis in [11] and determined important 
numerical parameters Ci and do which determines all the properties of the solutions. Furthermore, 
we have obtained expansion formulae for asymptotics at large distance. We have found that 
interesting two limits i) ^ m^, ii) m,^ ^ m^. The original linear sigma model reduces 
^2N'^-i non-linear sigma model in the i) limit and we have found a sort of non-Abelian global 
string solution there. In the second limit ii), we have obtained the SU{N) chiral Lagrangian 
coupled with a complex scalar field. We have also found a sort of Abelian global vortex solution. 
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It is a singular solution and can be identified with Abelian elobal vortex with an irrational U(l) 
winding number k = l/y N.\!] 

Here we give several discussions. The coupling between a global U(l) string and the U{1) 
Nambu-Goldstone boson can be constructed by using a duality between a boson and a two- 
form field Bfj^iy [2T]. In the same way the coupling of a global U{N) string and the U{N) Nambu- 
Goldstone bosons will be possible by using non- Abelian two- form |22j . 

In the presence of the U{1)a axial anomaly a term Vi($, <l>^) = c(det $ + det $''^) is induced, 
which gives a sine-Gordon potential to the phase. Then a U{N) vortex becomes a boundary of 
a domain wall [TOl [9]. In the presence of quark masses, a term V2($, = Tr[if($ + $"'')] exists 
with H a mass matrix. It remains as a future problem to study detailed structure of vortex 
solutions in the presence of these terms because the authors in [TOl [9] assumed constant profiles. 

In Subsec. 13.3.21 we have simply sent the mass to infinity to obtain the SU{N) chiral 
Lagrangian coupled with a complex scalar field s. However quantum mechanically we should 
integrate out the massive fields. This procedure generally induces higher derivative terms for 
remaining massless fields. The Skyrme term is such a term of the fourth order [23]. In our case 
with massless field s we will obtain the effective Lagrangian of the form 

= N\d,s\^ + \s\''TTid,Ud''U^)+^-^TTi[U^d^U,U^dM?)-V{\s\') 

= N{d^af + a^Ti{d^Ud^U^) + ^Ti{pd^U, U^d^Uf) - V{a^) (4.1) 

where e is a parameter determined by an explicit calculation and V is the potential in Eq. fl3.6ip . 

Note that there is no fourth order term for s. Because of the relation = s*d^s\i^ + 

|spt/^9^f/, the possible fourth order term $^9,^$]^ reduces to the Skyrme-like term in the 

Lagrangian (14. ip . When s is fixed to the vacuum expectation value, the Lagrangian (14. ip reduces 

to the U{N) Skyrme model and admits the usual Skyrmion solution for U{x) = expi[F{r)a ■ 

r/|r|]j^ It is interesting to note that the abelian vortex s ~ t;e*^ does not interact with the 

Skyrmion while a non-Abelian vortex does. It remains as a future problem to study interaction, 

scattering or absorption of baryons(Skyrmions) by non-Abelian strings. 

^ We have also investigated the Large N limit, and found that there exists a regular vortex string solution 

which has the same form with the usual U{1) vortex Eq. (j3.19p with winding number fc = 1 but with replacement 

7710 m^. See appendix for details. 

* The Skyrme model admits a topologically unstable string solution [24] which may be related to the pion 

string [5]. 
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In this paper we have studied U{N) vortices. Local and semi-local vortices with different 
groups [U{1) X G]/Z„„, where G is arbitrary group with the center uq [25], have been studied 
recently. In this framework the U (N) vortex corresponds to the case of G = SU (N) with rio = N. 
Global version of these vortices are also possible, especially the case of G = SO is related to 
vortices in the B-phase of ^He superfluids. 

Before closing this paper let us compare our global U{N) vortices with other types of U (N) 
vortices in the related models: 1) semi-superfluid U{N) vortices in high density QCD and 2) 
local U{N) vortices. In these models, the group structure is completely the same with the global 
case in this paper. However the energetics/interactions of vortices and the (non-)normalizability 
of the zero modes are significantly different. 

1) In high density QCD it is expected that color superconductivity is realized. There, the 
color symmetry SU{N)c and the flavor symmetry SU{N)p (with = 3) as well as the baryon 
U{1)b symmetry are spontaneously broken down to the color-flavor locked symmetry SU{N)c+f 
apart from the discrete symmetries. The corresponding vacuum manifold [SU{N) x U{1)]/Zn ~ 
U{N) is the same with that of the global U{N) vortices. In this case the SU{N) subgroup 
of the vacuum manifold U{N) is gauged and therefore only one massless Nambu-Goldstone 
boson for the U{1)b exists. The U{N) vortices here are called semi-superfiuid vortices [26]. 
In the asymptotic behaviour of the scalar field of a U{N) vortex, $ ~ f diag(e*^, 1, ■ ■ ■ ,1) = 

e*^diag ^e*"'^^, e~*^, ■ ■ ■ ,e~*^j, the latter non-Abelian part can be eliminated by a gauge 
transformation U{r, 6) = diag ^e"*'^^^*-''-*, e*^^*^''\ ■ ■ ■ , e*^^*-^'*^ with an arbitrary function F{r) 
satisfying the boundary conditions F{r = 0) = and F(r — > oo) = ij^ 

2) Next let us to compare the global U (N) vortices studied in this paper with the local U{N) 

vortices [12]. In this case too the symmetry breaking pattern is the same but a crucial difference 

^ This transformation is wcU-dcfincd because of the triviahty of the first homotopy group: ni[SU{N)] = 0. 
This transformation brings <I> to $ ^ ve^^rlN- This property makes the orientational zero modes of CP^^^ to 
be normahzable and to become the moduh (coUective coordinates) of the vortex [27| . In other words, the isotropy 
groups Hg at the infinities in the presence of a global U (N) vortex were physically different and depend on 6 
(although they are isomorphic), whereas, in the SU{N) gauged case, the isotropy group at the infinities in the 
presence of a semi-superfluid U{N) vortex is physically equivalent for any 6. The U{1)b is global and so the 
energy of a vortex remains logarithmically divergent. The asymptotic interaction between two semi-superfluid 
U{N) vortices is also essentially the same with the one between U{1) global vortices because of the above property 
[27l[T4]. Therefore it gives the universal repulsion between separated vortices. 
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is that the f/(l) symmetry is also gauged in addition to the color SU{N), and therefore there 
remain no Nambu-Goldstone bosons. Namely the vacuum is the unique at the infinity even 
in the presence of a U{N) vortex because the U{1) symmetry is gauged . Nevertheless we 
cannot gauge transform the asymptotic behavior of the scalar fields from $ ~ vc^'^In to Iat 
because we can define no regular gauge transformation well-defined in the entire space because 
of non-triviality of the first homotopy group: 7ri[f/(l)] ^ 0. Unlike the global U{N) vortices 
or semi-superfiuid U{N) vortices, these local U{N) vortices have finite energy because of the 
gauged U{1). The CP^~^ zero modes are normalizable. At the critical (BPS) coupling with a 
particular relation between gauge and scalar couplings, there is no static force among multiple 
vortices, with allowing the multi-vortex moduli space [28l 129] (see [30] for the moduli spaces 
of local U{N) vortices on a cylinder and a torus). Static interactions exist between vortices at 
non-critical (non-BPS) couplings. The force between two U{N) vortices was shown to depend 
on both CP^~^ orientations and positions |31j . 

If we gauge all the symmetry, the vortices are those in quiver gauge theories [32]. In this case, 
the diagonal gauge symmetry remains unbroken, and the vortices do not have orientations in the 
internal space. The final possibility which was not studied so far is the case that only U{1) is 
gauged. 
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A Large N Limit 

In this appendix we will derive an asymptotic form of the vortex solution in the large N limit. 
Before implementing the large N limit, we have to know how the parameters in the Lagrangian 
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scale in A^. From the observation of loop corrections (perturbation series in Ai and A2) to two- 
body meson scattering amplitude in a single channel, one can set Ai ~ 0{N~'^) and A2 ~ 0{N^^) 
in order to have a tenable perturbative expansion [33]. /i^ ~ because of no flavor degeneracy 
in a single channel. As consequences, one finds ~ thus r ~ 0{1). 

Now we are ready to see what happens in Large limit. After taking A^ ^ 00 as keeping 
and finite, Eq. (13.281) can be solved by 

g{r) = 1, (A.l) 

and the other equation (I3.27P becomes 

r + ^-^-^/(f-l) =0. (A.2) 

This is the equation for k = 1 Abelian global vortex. Note that is shown up in the equation 
and solutions are independent of rrifp unlike the case of Abelian global vortex given in Eq. (I3.19p . 
The solutions themselves happen to be identical to those for = (r = 1) with any finite 
N. 
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